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abstract 

Many systems whose evolution in time is governed by Partial Differential 
Equations (PDEs) are linearized around a known equilibrium before Computer Aided 
Control Engineering (CACE) is considered. In this case there are infinitely many 
independent vibrational modes, and it is intuitively evident on physical grounds that 
infinitely many actuators would be needed in order to control all modes. 

A more precise, general formulation of this grave difficulty (the spillover 
problem) is due to A.V. Balakrishnan [Applied Functional Analysis, Springer, 1981, p. 
233]. Let the system’s state vector x be an element of a separable Hilbert space X 
whose dimension is not finite; let A be a closed linear operator with domain dense in 
H which is the infinitesimal generator of a strongly continuous semi-group of 
transition operators Tit) for non-negative times t; and let B denote a bounded linear 
operator acting on another separable Hilbert space U (the control space) with range in 
H. Now consider the control problem dx/dt = Ax + £u, with x(0) given. Then according 
to Balakrishnan this system is not exactly controllable if B is compact. 

A possible route to circumvention of this difficulty lies in leaving the PDE in 
its original nonlinear form, and adding the essentially finite-dimensional control 
action Bu prior to linearization. In many cases it can be shown that the nonlinearity 
couples the system’s modes in such a manner that only a finite-dimensional subset of 
the modes is functionally independent, with the remaining higher-order modes 
nonlinearly dependent upon them. Hence control of all modes can be achieved by 
controlling only finitely many modes. 

One possibly applicable technique is the Liapunov-Schmidt rigorous reduction of 
singular infinite-dimensional implicit function problems to finite-dimensional 
implicit function problems. Such a procedure was employed by Leon Lichtenstein in the 
1930’s to prove the existence of a solution of the Navier-Stokes equations for a 

sufficiently small time-interval 0 ^ t < c. 

Omitting details of Banach-space rigor, the formalities of this approach are as 
follows. Let £ be a Fredholm operator with pseudo-inverse X; then there exist 
idempotent projection operators T = 9 - X#, Q = $ <£X whose ranges are 

finite-dimensional and such that a NASC for l*) £x - 5(x) = 0 is that (**) x = Tx + 
X9Hx) 8c Q^(x) = 0. Thus one may set x = u + v where v satisfies the auxiliary equation 
v = X^(u + v) and u the bifurcation equations [Verzweigungsgleichungen] Qf(u + v) = 0. 
Typically one may solve the auxiliary equation by (contraction) iterations to find v — 
£( u ) where now v is infinite-dimensional but u is finite-dimensional and then insert 
the result into the (finitely many) bifurcation equations to define a 

finite-dimensional vector function flu) = Q^lu + &lu)) such that (*) is equivalent to 
/(u) = 0. In summary, a NASC for (*) is 

(***) x = u + £(u), £u = 0, PSlu) = 0, Tu = u, flu) = 0. 

As an illustration the auxiliary equation and bifurcation equations for the 
problem of deflection of an in tension (a Q > 0) EXTENSIBLE beam (a j > 0) is considered, 

including viscous damping (a^ > 0) and Balakrishnan-Taylor damping (a^ > 0). Here 


u + a *u = a *u + 

tt 2 xxxx 3 xxt 


( r 1 " 2 

r pL -i 2(n+6)+l> 

a + a • (u ) -dx + a • 

l 0 1 J 0 x 4 



* U 


V t £ 0, 0 ^ x ^ L. As the dimension N of the bifurcation equations increases, the 
result approaches an W-dimensional truncated eigenexpansion (provided that the initial 
deflections and their initial spatial and temporal rates of change are not too large). 
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Pref ace 


The basic idea behind the present paper is simply; 


SUGGEST 8 GIN 

Don’t linearize a PDE until after its reduction to a finite-dimensional ODE. 


This idea can be implemented by means of the following analytical procedure: 


L D APU1N0 V - S OHM 0 DT BIFURCATION EQUATIONS: 

A rigorous reduction of a singular inf in i t e -d i me ns i on a 1 implicit equation to 
the problem of an equivalent, merely / in i t e -d imens i onal i mplicit equation. 


This suggestion is presented as a possible technique for circumvention of the 
famous "Spillover Problem." 


Introduction 


If the problem of control of a flexible structure is linearized before one 
considers the control aspects, then frequently it leads to an abstract problem in 
functional analysis of the type of the following system of ordinary differential 
equations: 

~ = Ax •¥ Bu t x(0) = x°. (1) 

Here x e H is an element of an infinite-dimensional state space taken to be a 
separable Hilbert space. Also u € 1/ is an element of the control space, taken to be 

another separable Hilbert space. We take A: D -» K°° to denote a closed linear mapping 

of the dense linear-subspace domain D into H°° = X> which is the infinitesimal generator 
of a strongly continuous semi-group of transition operators T(t) for t > 0. Finally we 

require that B: 1/ H be a bounded linear operator. 

The celebrated "Spillover Problem" now has an exact formulation by means of: 

THEOREM. (Balakrishnan, [5], p. 233.) If B is compact, then if is NOT exactly 
controllable. 

A compact operator is one which can be approximated arbitrarily closely by an 
operator whose range is finite-dimensional. Therefore the practical import of the 
preceding theorem can be phrased as: if a linear system has infinitely many 
independent modes of motion , it cannot be controlled completely with a finite- 
dimensional actuator suite . 

This suggests that complete control of a flexible structure by a finite actuator 
suite is foredoomed to impossibility. However, there may be a way to circumvent this 
difficulty. Note that the preceding theorem has been proved only in the case that the 
dynamical system if is linear. The purpose of the present paper is to demonstrate that 
for many flexible structure problems, such nonlinear mechanisms as Balakrishnan-Taylor 
damping will couple the higher order modes of motion to the lower order modes in such 
a way that only a finite number of the lower order modes is functionally independent. 
This suggests that a finite actuator suite could control such a system. However, we 
defer consideration of the control problem and deal here only with the free motion of 
an uncontrolled, but intrinsically nonlinear system. Our purpose is to stimulate 
further research into this approach rather than to present a finished theory. 
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Liapunov-Schmidt Bifurcation Theory 


For the reader’s convenience we recall the salient features of this theory from a 
purely formal point of view. The details of Banach space rigor can be found in pages 
173-177 of Deimling [2] and other texts on nonlinear functional analysis [3], [4], 

Let £ denote a Fredholm operator, which may be singular, i.e. there may exist 
elements u * 0 such that £u = 0. Let K denote a pseudo- inverse of £, i.e. a linear 
operator such that 


KUK 

= K , 

i eK£ = 1. 

(2) 

Now define projection operators 

? = 5 

- X£, 

Q = £ - £K\ 

(3) 


it is readily verified that P Z = P, Q 2 = Q, i.e. these operators are idempotent, 
which justifies referring to them as projection operators. Note that P is a right zero 
of £, and Q is a left zero of £. 

Let (F = 3H - } denote a nonlinear operator. Then it is easy to verify the 


equivalence of the implicit equation problem 

(*) £x = ?{x}, £P = 0, Q£ = 0, (4) 

and the problem 

(*«) x = Px + X?{x>, Q9Tx} = 0. (5) 

Now define u = Px and verify that Pu = u; then we can replace (**) by 

(definition) X = U + V, (6) 

(AUXILIARY EQUATION) v = K&{u + v>, (7) 

(BIFURCATION EQUATION) Q2{u + v> = 0. (8) 

Another name for the Bifurcation Equations is Branching Equations. 


Suppose that the right-hand side of (7), regarded as a function of v, has a 
global Lipschitz constant less than unity. Then by the well-known principle of 
geometric convergence of Contraction Mappings we may, for each fixed u, define a 


nonlinear mapping v = £?{u} as 

v = lim v k , v k+1 = + v k >, v° = 0, ( k = 0, 1, 2, 3, ■ • • ). (9) 

k->oo 

Here £* is the resolvent of the auxiliary equation in the sense that 

£ = Knu + £}. (10) 

Hence we may eliminate the auxiliary equation and replace v in the bifurcation 
equation by S* to obtain a new finite-dimensional equation 

fiu) = Q?{u + £{u» = 0, (ID 

which is equivalent to the original infinite-dimensional implicit equation. Thus 
(*) <^> (**) <=> (***) x = u + ^{u}, £u = 0, T^{u) = 0, Tu = u, f(u) = 0. (12) 

If the original functional equation was analytic, then the final 

finite-dimensional equation f(u) = 0 will also be analytic. 


If £ was non-singular, then K = £ \ whence T = 0, Q = 0, u = 0, the bifurcation 
equation does not arise, and the resolvability of the auxiliary equation is equivalent 
to the resolvability of the original equation: 

£x = 5{x) <=* x = £ l ^{x) = (13) 

Finally, if was singular, then the linear part F = f (0) = (df /dx ) of fix) is 

X I J 

necessarily also singular. Typically then the solutions of fix) = 0 will not be unique 
and one studies the branching of these solutions by such methods as Newton's polygon. 
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Deflection of an Extensible, Nonlinearly Damped Beam 


Let u denote the normal deflection from equilibrium. Then the vibrations of the 
beam can be described by u = u(t,x), 0 ^ x < L, 0 ± t < +co, which satisfies the PDE 


p*u + £J*u = C-u + I H + 

tt xxxx xxt I L 


( H * T 5 '* * r * ) 


(14a) 


4 > - ( 1 / 2 ) 


■ f {u^) Z -dx, 


(14b) 


^ (u U ^)*d>:j 


2(n+B)+l 


( 0 S 0 < 1/2 ), ( n = 0, 1, 2, 3, 


where 


t € 1R = [ 0 , +00 ) = { T I 0 ^ T < +oo }, 

+ 1 

x <= J = [0,1] = { £ | 0 < e s 1 }, 
p = density, 

E = Young’s modulus of elasticity, 

I = cross-sectional moment of inertia, 

C = coefficient of viscous damping, 

H = axial force (tension or compression) 

Ac = cross-sectional area, 

L = length, 

T = Balakrishnan-Taylor damping coefficient. 

Now define the constants 

a^ = H/p, a^ = EAc/{2pL) t a - EI/p t a ^ = C/p, a^ = T/p 
Then the PDE (14) can be expressed as 


u + a • u = a * u +|a 

tt 2 xxxx 3 xxt 


(14c) 


(15) 


~L r~L -i 2{n+B)+l\ 

+ a • (u ) Z -dx + a •II (u u )-dx -u . (16) 

0 1 J x 4 Un XXt j Jxx 


Boundary Conditions 


As usual, we require that 

u(t,0) = 0, 


u(t,L ) = 0, 


(17a) 


u (t,0) = 0, 

XX 


u (t,L) = 0. 

XX 


(17b) 


Initial Conditions 


Let <f> and \p be functions of x defined on J with the following smoothness requirements. 
The function <p should be continuously once differentiable on J and its second 


derivative <p” should exist almost everywhere on J and be [Lebesgue] square-integrable 
on J. The function should be continuous on J, which of course implies that it is 
[Riemann] square-integrable on J. These smoothness requirements may be summarized as: 

„<i) 


<t> € 
\p e 


<P 9 


C ( 0 ) (J), 


0 € 


LU), 

l 2 ui 


(17c) 


(17d) 
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(17e) 


Now 0 and 0 are used to define the initial conditions on u as follows: 

u(0,x) = <p(x), u^( 0,x) s 0(x), 

0(0) = 0(L) = 0, 0(0) = 0U) = 0. 


(17f ) 


Normalized Constants 


For future convenience, we define 


b = (n/L) 2 -a , b = (ti/L) 4 - (L/2)-a , b = (re/Z.) 4 -a , 
0 0 1 12 2 

, /n 2 , , /r v4U+n+fi] , ,2(n+B)+l 

b = (n/L) -a , b = (tt/L) -(L/2) -a . 

3 3 4 4 


Function Space Coordinatization 


(18a) 

(18b) 


Define the complete orthonormal set on L^fJ) by 


e - e (x) = sin{kn[x/L]), ( k = 1, 2, 3, 

k k 


). (19) 

Assumptions (17c, d) imply that there exist (as l.i.m.) sequences {a^}, such that 

00 

(20a) 
(20b) 


0 = 0(x) = ^a^-e^tx), 
k = l 
00 

ip = \p{x) = y P *e (x), 
U k k 


E *S 2 

= 1 

00 

l *1 


< + 00 , 


k = l 
oo 


< +oo. 


k = l k = 1 

Now we can seek to find a sequence of time-varying functions {u (£)} such that 


u(£,x) = ) u (t) e (x), 

L k k 


k = 1 

where, by (17f), the initial values and initial rates of the {u (£)} must satisfy 

k 

u (0) = a , u (0) = /3 , ( k = 1, 2, 3, • • • ). 

k k k k 


( 21 ) 


( 22 ) 


Infinite System of ODEs 

Insert the series expansion (21) into the PDE (16) and use the orthonormality 
property of the complete basis {e^} to derive an equivalent infinite system of ODEs: 

ii + b k 2 -u + (b + b $ + b 0 + b k 2 )-k 2 -u =0, ( k = 1, 2, 3, • • • ), (23a) 

k 3 k 0 1 4 2 k 


1 ■ l ■ 

• ■ [I 'Vj 


2(n+fi)+l 


(23b) 


(23c) 


where the solutions of (23) are required to satisfy the initial conditions (22). As an 
alternative to (23a), in which the linear and nonlinear terms are displayed 
separately, we may write 

ii + b k 2 -u + (b + b k 2 )-k 2 'U = -$*k 2 *u , ( k = 1, 2, 3, * * * ), (23d) 

k 3 k 0 2 k k 
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$ 


b $ + b 0 

1 4 


oo 




.2 

J 


U -u 

J JJ 


2(n+fl)+l 


Energy Integral 


(23e) 


Multiply each equation of (23a) by 2u^ and sum over all k t using the identities 
2£i *ii = (u ) and 2u *u s (u ) , to obtain 

k k k k k k 

y u 2 + 5 * y k 2 *u 2 + b • y ^ 4 - u 2 + (b /2)- f y k 2 -u 1 + 

A k 0 k=i k 2 a k 1 Lfc k J 

+ 26 3 J‘ [jV-i 2 ] dt ♦ l V 2 i f [i * 2 -u -u ] 

J 0 L k = l J J 0 L k = l J 


i 2(n+fi)+l 


dt 


= & = constant = 


K 


(24a) 

(24b) 

(24c) 


■ 4 -f * 


k =1 


2 u 2 
dx + b •— j- 
o L 


y (b + b k 2 )-ic 2 *a 2 + (b / 2 )|" y k 2 -a z 1 = 

L 0 2 k 1 L kj 

k — 1 k — 1 

L L L 2 

f (<p’) 2 dx + b 2 -ff (0”) Z dx + (b/2)[?-f (0’) 2 dx] s (24d) 

J 0 ^ 0 1 L J o J 


= Ifr 2 + b *(0’) 2 + b *(0”) 2 +■ (b 

0 2 


- -» 2 

/ 2 )- ( 0’) 2 I , 


(24e) 


where in the last expression we have used an obvious notation for the mean-square 
values of 0, <f>\ and <p" on J = [0 t L]. 

A key technique in what follows is the use of (24) to obtain a priori bounds on 
the solutions of (23). 


Vector Notation 

For convenience we shall denote the infinite column whose rows are ku by x, and 

k 

we shall partition x into a finite-dimensional component u and an infinite-dimensional 
component v as follows: 

X = ^ v ] ^ ( ^ U k ^ u e (R N , v € (R 00 , (25a) 

u = ku , v = kv t { v = u , k £ N ). (25b) 

k k k k k k 

A Priori Bounds 

Obviously 


00 


iixii = null 2 + 

while from (24a) it is easy to infer that 

„V„ 2 S j> 2 - u 2 , 

k = l 

(26) 

(b + b )■ llxll 2 + 
0 2 

(b/2)- llxll 4 < 8 , 

1 0 

(27) 


which implies that 
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Ilxtl < R, (28) 

where (choosing the numerically stable quadratic root formula) 

R = (2& Q /{(b o + b 2 ) + [(b Q + b z ) 2 + Zb& Q ] l/2 )) l/Z * (g/Ib, + b 2 ])1/2 ' (29) 

Similarly, 

(W+l) 2 - II v II 2 3 (N+l) 2 - £ k 2 -v 2 S £ k 4 -v 2 3 - I k4 ‘ u k " e o /b 2 ' (30) 

k-N+1 k = N + 1 k=N + 1 k = l 

and 


Yu 2 * & , 
L k 0 


which together with (30) yields, via the Cauchy-Schwarz inequality, 
00 , , QO \ ✓ oo \\ 1/2 . „ 


00 r f w \ f w \ 

1 * { H • I‘ 4 '\ 

j=l v ^ k = l 7 V k = 1 7 


s (g [g /b ]) 1/2 = g /(b ) 1/Z . (32) 

0 0 2 0 2 


Finally, using (29) and (32), 


J = b • llxll 2 + b • (x-x) 2(n+fi)+1 ^ b *R 2 + b • [£ /(b ) 1/2 ] 2(n B) < b •© , (33) 

1 4 1402 SO 


b s [b /(b + b )] + b /(b ) 

5 1 0 2 4 2 


n+fi+(l/2) 


where we have used the assumption that 

& < 1. (35) 

o 

Now let B N (p) and B°°(p) denote balls of radius p in ER and (R°° where x = u®v is 
considered to be an element of the Cartesian product IR ®1R ; we shall show that for 
all sufficiently large values of N, B N (R^)®B (R^) is a subset of the ball defined in 
tR N <s>(R°° by (28), where 0 < c < 1 is an arbitrarily small positive constant, where 


and where 


In fact, 


R = R • (1 - e) , 
e 

R = (g /b ) 1/2 /(W + 1), 

N 0 2 


N + l * 18 /b ) l/2 /(Rc 1/2 ). 

0 2 


Hull 2 + llvll 2 S (i? ) 2 + (R ) z S g 2 -(l - e) + R 2 -e = R z . 
e n 

Initial Conditions 


Define 


and note that 


Next, define 


a = ( k-a ) € 1R N , a (N) = ( [N + k] *a ) e OR , 

k N+k 


a = u(0), a = v(0). 


b 3 ( p ) 6 R n , b <N) 3 ( 0 ) e r", 

k N+k 
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and note that 




b (N) = © 1 • v(O). 

N 


b = © -u(0), 

where the bounded linear operators D~ l and © _1 are represented by matrices 


© 


= diag ( 1/Jc ), ©^ = diag ( l/[k + N] ). 


(40d) 


(40e) 


Resolution of Linear Part 


u + b k 2 -u + (b +b k 2 )-k 2 -u = - 
k 3 k o 2 k 


Alternatively, in vector notation 


So 1 * + b 3 -I) o X + (b o + V®o ),JJ o x = " { b r 11x1,2 + b 4 -(x-x) 
b Q = axial force coefficient, 




,2(n+B)+l 


3, 

), as 

2 

: • u . 
k 

(41) 

) X, 
0 

(42) 


b = extensibility coefficient, 
b = elasticity coefficient, 
b - viscous damping coefficient, 


(43a) 
(43b) 
(43c) 
(43d) 

= Balakrishnan-Taylor nonlinear damping coefficient, (43e) 

where the unbounded linear operator © q is defined by the infinite matrix obtained by 
taking IV = 0 in the reciprocal of (40e). Upon multiplying (42) through by ©^x and 
integrating, we obtain the energy integral (24) in the form 

i: 


ll2) o x!l + b o‘ 11x11 + (bj/2) • II x II 4 + b^ II Ex II 2 + 2 I {b^llxll 2 + b • (x • x) 2(n+fi+1) > dt = 

z o 4 


= & = ilb°°l| 2 + b ’ lla°°ll 2 + (b /2) • fl a°°ll 4 + b *112) a°°ll 2 , 

0 0 1 2 0 


where as before 


00 N ,00 __ , N 

a = a®a , b = b®b . 


For future use define c and d bv 

N N J 


IN + l) 2 - lia N li 2 


HbV = 


uu 

E 4 2 2 

k • a = c -> 0, N -» +oo, 

k N 

k =N + 1 
oo 

l 0* 5 S l * °. * -> 


(44a) 

(44b) 

(45a) 

(45b) 


k = N+1 

Homogeneous Linear Part 


Henceforth we shall assume that the system is underdamped, i.e. that 

b < 2(b )’ /2 

3 2 

The characteristic polynomial of the homogeneous part of (41) is 


(46) 
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A 2 + b k 2 -X + (b +b k 2 

) - k 2 = 0, 

(47a) 


k 3 k 0 2 



which has roots 

X = -p ± i * V , 

( i 2 = -1 ). 

(47b) 


k k k 




= k 2 -p Q , ^ = 

X. 

3 

CM 

II 

(47c) 

hi 

o 

b /2 , it) s (i/2) ■ {4b - 

3 k 2 

- (b ) 2 + 4(b /k 2 )> 1/2 , 

3 0 

(47d) 

where, obviously, 

ji / 0) £ b /{4b - 

(b ) 2 } l/Z . 

(47e) 


0 k 3 2 

3 


From (47a) the general solution of the homogeneous linear part of (41) is 

of the form 


u = exp(-u t) • {A - 005(12 t ) + B -sin(v £)}, 

k k k k k k 

(47f) 

where A and B are arbitrary constants. Specializing these constants 

Lr \s 

in order to use 

Lagrange * s var i at i on 

of constants formula to 

re-express (41) in terms of an impulse 

response convolution 

with the right-hand side we 

get 


u (t) = 

k 

d (U (t)-a + (d l2) (t)/k) •£ - 
k k k k 

f d <3) (t - x) $ li(x) dx, 
J k k 

J n 

(47g) 


where 


d (1) (t) = exp(-/i t)'{cos(v t) + (u /u )'sin(v £)}, 

l( k k Ok k 

d i2 \t) s exp(-p t)-{[l/Uc-u )]-sin(i> i)>, 
k k k k 


d (3> Ct) = k-d (2 \t). 

k k 

Next, define three finite and three infinite diagonal matrices as 

D = diagi d*^) : [R N -> (R N , ( J = 1, 2, 3 ), ( k = 1, 2, 3, * • • , N ), 

J k 

D = diag( d (J) ) : (R°° -» R°° , ( j = 1, 2, 3 ), ( k = 1, 2, 3, • ■ • ). 

j,N k+N 

Finally, multiply (47g) through by k in order to convert to vector notation: 

x(t) = D (t)-a“ + D (t)-b“ - I D (t - t)-»-x(t) dr, 

h * ^ o 

- 

x(t) = D ( £ ) • a°° + D ( t ) - b°° - i> (t - x)- 4 >-x(x) dx, 

1,0 2,0 J 3,0 


(48a) 

(48b) 

(48c) 

(49a) 

(49b) 

(50a) 

(50b) 


<J> = b 


b -(x-x) 

4 


2<n+G)+l 


(50c) 


(In deriving (50b) we used the fact that q ( 0) = 0.) The fact that we have been 

able to reformulate the original PDE boundary-value and initial-value problem in the 
form (50) is the equivalent of (13), i.e the non-singular case, wherein there is no 
requirement for a Bifurcation Equation since the entire system now has the form of the 
Auxiliary Equation. If now we can prove that the iteration, for m = 0, 1, 2, 3, ■ , 
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(51a) 


x m+1 (t) = D (t) • a°° + D (t) • b° 

1,0 2,0 


J £> 3 Q (t - T)-*(x" 1 ,x m )-x m (T) dr. 


x°(t) = 0, x J (t) H D (t)-a°° + D (51b) 

1,0 2,0 v ; 

converges, then we have constructed a solution of the original problem. 

To avoid certain difficulties, we shall consider this iterative solution only for 
the problem 

b 4 = 0, = $(x) = b^llxll 2 ; (52) 

the more general problem will be approached by a non-constructive homotopy method. 
Consider now the fixed-point problem 


x = 3Tx}, II x II = sup llx(t) II , 

CO 

(t>0) 

? = 3 c {x(*)}(t) = D iQ (t)-a + - J Q (t - t)*$(x(t))-x(t) dr. 

We want to find a Lipschitz constant for 3% i.e. a constant x such that 

H3{x J+1 > - ^<x J }ll oo < x-llx J+I - x J ll m V x J e B°“(R), ( j = 0, 1, 2, ••• ). 
Later we shall prove that 


f + “ 

IID (x)ll dx < b = 4/{b (4b -(b) 2 ) 172 }. 

J 3,0 6 3 2 3 

In the Hilbert space norm of (R°° it is clear that 

II { llx II } x - (llx II } x II — II {llx 2 ll 2 } (x 2 - x 1 ) + {(x 2 + x 1 )-(x 2 

s {llx 1 !! 2 + llx 2 |l 2 + llx 1 II • llx 2 H> - llx 2 - x 1 !! s 3J? 2 -llx 2 - 
Hence we may take 


- x 1 )} x 1 
i„ 


(53a) 

(53b) 


(54) 

(55) 


(56) 


k = 3R Z b^b b i 3{b i -b & /[b o + bj}^ - & < 1, (57a) 

if & q is taken to be sufficiently small that 

e < (1/12) -0-<(b + b )-b -[4b - (b ) 2 ] 1/2 }/b . (57b) 

^ 0 2 3 2 3 1 

The a priori bounds (28M29) apply to the first iterate (51b) and so the first 

iterate is inside the ball B°°(R); now, using (54) with j = 0, 1, 2, etc. it is clear 

that x 2 , x 3 , all remain in the ball provided that llx'll + tf-llx 1 !! + 0 2 -||x 1 ll + ■■• < R 

00 00 00 ’ 

i.e. summing the geometric series, provided that 

^ Hx 1 Hoc " 0)> < R, (58a) 

which is readily obtainable simply by taking ■& so small that 

& < 1 - (llx 1 ll (n /R). (58b) 

Now we can apply the well known principle of contraction mappings (also called the 
Banach fixed point theorem and Caccioppoli's fixed point theorem) to prove that the 

map (53) has a unique fixed point in the ball B “(K) and that this fixed point can be 

computed constructively by the iterative procedure just described. 
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Next, replace the coefficient b^ by pb^, where 0 * p s 1, and note that the a 

priori bounds previously derived remain valid for all values of p € [0,1]. Hence we 
may infer the existence of a solution for all p € [0,1] by a homotopy method. For 

sufficiently small values of p, this is Poincare's method of analytic continuation of 
solutions of functional equations. Here the fact that the problem was non-singular at 

p = 0 implies that the Leray-Schauder Index [2], [3] or topological degree of the map 

is of magnitude unity at p = 0. Consequently the existence of the a priori bound for 
any solution on 0 ^ p ^ 1 implies that there is a continuum of solutions connecting 

the solution at p = 0 with one at p = 1. Thus in summary we have proved the following 
result. 


Existence Theorem 


THEOREM . If </>\ <)>” > and \p [as defined in (24d-e)] are all sufficiently small, 

then the nonlinear functional PDE (16) has at least one solution which exists for all 
t > 0 and satisfies the boundary conditions (17a-b) and initial conditions (17e) as 
well as the a priori bounds 


2 

L 


l 


[u (t,x)] 

X 


2 


dx 


R 


8 /(b + b ), 
0 0 2 


(59) 


where 



[u ( t,x )) 2 



dt 


s S /(2b ), 
0 3 


8 

o 


ip 2 + b Q -(<p') 2 


+ b •( ip ”) 2 + (b /2) • 
2 1 



(60) 


(62) 


Proof . Existence has already been proved. Equation (59) holds because 

the left-hand side of (59) is by (26) equal to llxll 2 which in (28) is proved smaller 

than Ft 2 . Equation (60) follows from inspection of (24a) and (24b), which hold for all 
t a O, making it permissible to let t -» +». 


Rigorous Truncation 


In "naive truncation" one simply sets v = 0, i.e. in (41) one takes 

u (t) s 0, ( k = N + 1, JV + 2, N + 3, • • • ). (63) 

k 

Here we shall prove that the solution proved above to exist can be derived as in the 
Liapunov-Schmidt bifurcation theory, and that, for all t ^ 0, 

II v( t ) II ^ R *exp(-b [N + l] 2 -t/2) ^ Rc 1/2 t (54) 

N 3 

provided that N is taken larger than the lower bounds in (38) and in (68b)-(69) below. 

By inspection of (50), we can express the problem as follows. 


Bifurcation Equation 


u(t) 


D (t) ■ a + D (t) - b 
1 2 


f D (t - x)-l-u(x) dx, 

^ o 3 


(65a) 
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$ = b ■ {Hull + llvll 2 } + b *{u*u + vv) 2(n+flM 

i 4 

Auxiliary Equation 


(65b) 


v(t) 


D (t) • 

1.N 


D (£)■ 

2»N 


- f D 

L 3 > N 


U - x)-$*v(x) dx. 


More A Priori Bounds 


(65c) 


From (48a), 

Id^’l s 0 + (M 0 /w k ) 2 } 1/2 -exp(-/i o A 2 t) s {1 + i> 2 /[4b 2 - (b 3 ) 2 ]} 172 -exp(-b 3 i/2) s 


= [1/{1 - (b 2 / 4b )} 1/Z ] -exp(-b t/Z). (66a) 

3 2 3 

The induced Euclidean norm of any diagonal matrix is equal to the absolute value of 
its (absolutely) largest element. Hence 

H£HI ^ [1/(1 - (b 2 /4b 2 )} 1/2 ] *exp(-b t/2). (66b) 

Similarly, 

IID 2" " iSgifn K 2) ! " {l/(*u k )>-exp(-b 3 t/2) * 

- [2/{4b z - (b 2 )} 1/2 ]-exp(-b 3 t/2), (66c) 

" D 3 N " < 1/ (^ k »-exp(-b 3 t/2) s [2/{4b z - (b 2 )} 172 ] • expi-b^t/2), (66d) 

HD N « * 0/(1 - (b 2 /4b z )] 1/2 )-ejcp(-b 3 [W + l] 2 t/2), (66e) 

IID 2,n" - <[2/{4b 2 - ( b 3 )) 1/2 ]AN + l)>-ejcp(-b 3 [W + l] Z t/2). (66f) 

HD 3N II s [2/{4b z - (b 2 )) 172 ] -ejcp(-b 3 [W + l] 2 t/2). (66g) 

Integrating (66d), one obtains 


Similarly, 



II D (x ) II dx 
3 


< b 

6 


4/{b (4b - (b ) 2 ] 172 }. 

3 2 3 


(66h) 


r + “ 

HD (t)II dx £ b AN + l) 2 , 

J Q 3,N 6 

which we have used above in (55) in the case N = 0. 

Next, from (65c), (66e,f,g), and (45a, b), 

Hv(t) 11 s y -exp(-b IN + l] Z t/2) + (b b S /(N + 1) 2 >- II v(t) II, 

where 


(66i) 


(67a) 


y = (1/(1 - (b 2 /4b )] 172 >lla N ll + {[2/{4b - (b Z )> 172 ]/(W + l)>llb N ll < 

r* 3 2 2 3 
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£ |{1/[1 - (b Z /4b )] 1/2 }-e + [2/{4b - (b 2 )> 1/2 ]-5 \/(N + 1) s cr AN + 1). (67b) 

[ 3 2 N 2 3 N J N 

Consequently, if is so small that 

b b & < 1/2, (68a) 

6 5 0 

i.e. so small that 

S < l/f<8/(b [4b - (b ) 2 ] 1/2 }-{[b Ab + b )] + b /(b ) n+B+(1/2) >] (68b) 

0 [323 102 42 J 

then we may subtract the second term on the right hand side of (67a) from the left 

hand side, and then, by (45a, b), making N still larger (if necessary) so that 

O' /{I - (1/2))} 3 20- £ (S /b ) 1/Z 3 (N + 1)-R (69) 

N N 0 2 N 

we can ensure that, for all t i 0, 

7 i /2 

llv(t)ll £ R -ex p(-b IN + 1] -t/2) £ R £ R e , (70) 

N 3 N 

as claimed in (64). 

Quite similarly, from (65a) and (68a), 


Hull < II D II * Hall + HD II • lib II + (b b & )• Hull < HD II -Hall + ll£> II • llbll + (1/2) • Hull , (71a) 

1 2 6 5 0 1 2 

whence, precisely as before, 

{1 - (1/2)}- Hull ^ It D II * Hall + IID IMIbll, (71b) 

i.e. 

Hull ^ 2* {IID II * Hall + IID IMIbll} < R , (71c) 

1 2 c 

provided only that, by (66b, c) 

Hall £ (1/4) -{1 - (b Z /4b )} 1/Z -R , (72a) 

3 2 G 

llbll £ (1/8) • {4b - (b Z )> 1/Z -R , (72a) 

2 3 G 

which can be assured simply by taking <p y and 0 small enough. 

Recalling that, by (62), & can be made arbitrarily small by making the initial 

mean-square spatial rates of change (p\ 0”, and ip of u , and arbitrarily 

small, we may conclude that when those initial rates are sufficiently small then every 
solution of (65a, b) must, for all t ^ 0, satisfy a priori 

1 /9 

II u( t ) It < R = R’ (1 - e) , (73a) 

G 

and every solution of (65b, c) must, for all £ £ O, satisfy a priori 

II vtt) II £ R £ R-e 1/Z , (73b) 

N 

whence finally, as in (39), every solution of the combined systems (65a, b,c) must for 
all t £ 0, satisfy a priori 

II x II 5 <llu(t)li z + llv(t)ll Z } 1/Z £ R. (73c) 

Conclusion: Naive vs Rigorous Truncation 

Let u(t), u(£) be arbitrarily given continuous N-vector functions and insert them 
into the nonlinear function $ defined in (65b), and then insert this functional of v 
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only into (65c), giving an infinite system of integral equations for v(t). As in (28) 
through (35), the nonlinear term can be made arbitrarily small by taking 

arbitrarily small. Similarly the nonlinear term can be made to have a Lipschitz 
constant less than unity by restrictions upon as in (53)-(58). 

This infinite system of Auxiliary Equations can be solved by one of the methods 
illustrated above (iteration or homotopy), and the result inserted into the 

Bifurcation Equations (65a) to provide a finite-dimensional system of functional 
integral equations exactly-equivalent to the the original infinite-dimensional system. 

When the arbitrarily given u(t), u(t) in the Auxiliary Equations (65b, c) are 

taken to be the projection into [R N of the solution x(t) e (R N «^°° proved to exist in the 
Theorem concerning (59)-(62), then the Auxiliary Equations have a solution 

corresponding to the projection v(t) of x(t) into R°°. The resulting Bifurcation 
Equations then must be satisfied by the same u(t) used to define the Auxiliary 
Equations. However, all of the a priori bounds proved above to apply to the solution x 
of the complete problem now apply to the projections u, v of the rigorously truncated 
problem. 

Consequently we can compare the rigorous version of the Bifurcation Equations, 
namely (65a), with the naively truncated version wherein one sets v = 0, and note that 
that, as becomes sufficiently small for the bounds to apply, then as N becomes 

arbitrarily large the difference between the solutions of the naively truncated 
version of (65a) and its rigorously truncated version becomes arbitrarily small. 

Consequently in attempting to solve the given nonlinear function PDE boundary- 
value initial-value problem, we may be confident that if we truncate naively for some 
finite N , the results become arbitrarily accurate as N increases without limit 

provided that the initial mean-square spatial rates of change <p * , <£”, and \ft of u , 

X 

and are kept sufficiently small. 

Further research is needed in order to ascertain whether or not this conclusion 
would still apply if the homogeneous boundary conditions were replaced by 
inhomogeneous boundary conditions corresponding to a finite number of actuators. 
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